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Abstract. We have calculated the evolution of low metallicity red giant stars under the 
assumption of deep mixing between the convective envelope and the hydrogen burning 
shell. We find that the extent of the observed abundance anomalies, and in particular 
the universal 0-Na anticorrelation, can be totally explained by mixing which does not 
lead to significant helium enrichment of the envelope. On the other hand, models with 
extremely deep mixing and strong helium enrichment predict anomalies of sodium and 
oxygen, which are much larger than the observed ones. This latter result depends solely 
on the nucleosynthesis inside the hydrogen burning shell, but not on the details of the 
mixing descriptions. These, however, influence the evolution of surface abundances with 
brightness, which we compare with the limited observational material available. Our 
models allow, nevertheless, to infer details on the depth and speed of the mixing process 
in several clusters. Models with strong helium enrichment evolve to high luminosities 
and show an increased mass loss. However, under peculiar assumptions, red giants reach 
very high luminosities even without extreme helium mixing. Due to the consequently 
increased mass loss, such models could be candidates for blue horizontal branch stars, 
and, at the same time, would be consistent with the observed abundance anomalies. 
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1. Introduction 



The observed anomalies in CNO-, NeNa- and MgAl- elements in globular cluster red 
giants (see Kraft 1994 and Da Costa 1998 for reviews) are unexplained in canonical low- 
mass star evolution theory and indicate effects beyond the standard picture. At least 
for anomalies in those isotopes participating in the CNO-cycle models relying on the 
assumption of an additional, non-standard mixing process inside the stars have been 
presented, which explain convincingly the observatio ns, including the e volution of the 



carbon abundance along the RGB, i.e. with time (see Charbonnel 1995; Denissenkov fc 



Weiss 1996 ; Cavallo et al. 1998 ). This mixing is supposed to set in after the hydrogen 



burning shell has reached the composition discontinuity left behind by the firs t dredge- 
up on the red giant branch (RGB), i.e., after the so-calle d RGB bump (e.g. Sweigart 
fc Mengel 1979 ; Charbonnel 1995| ; [Charbonnel ct al. 1998 ) where the molecular weight 
barrier between hydrogen burning shell and envelope is at a minimum. It is usually 
described in terms of a diffusion process of certain efficiency and penetration depth. 
The CNO anomalies and their correlations with brightness and metallicity (which are 
expected qualitatively from nucleosynthesis arguments in standard models, as discussed 
by Cavallo et al. 1998) are thus explained by a purely evolutionary picture ( ^mith fc Tout 



1992; Denissenkov & Weiss 1996; Boothroyd & Sackmann 1999a). The physical origin 



of the mixing process is believed to be found in differential rotation of the star and the 
parameters used in some of the presently available calculations have been derived from 
existing theories (e.g. Zahn 1992; Maeder fc Zahn 1998| ), which are, however, far from 
being complete. 

A similar situation holds for oxygen and sodium, which are found to be anti-correlated 



in globular cluster red giants (Kraft et al. 1993). Denissenkov fc Denissenkova (1990) 



and Langer et al. (1993) showed that this could result from the mixing of elements 
participating in the ONeNa-cycle, which operates at higher temperatures than the CNO- 
cycle and therefore requires deeper mixing. Denissenkov fc Weiss (1996) demonstrated 
how all anomalies of the mentioned elements known at that time can be explained by the 
deep mixing scenario. Their calculations were done by using canonical red giant models, 
which evolved along the RGB without any mixing, and performing the mixing and nuclear 
reactions in a post-processing way. For this approach to be correct it is necessary that 
the evolution of the background models is not affected by the mixing process. In fact, the 
mixing necessary to reproduce the observed anomalies was always so shallow that only 
very small amounts of hydrogen/helium were mixed between envelope and shell, even 
in the case of the 0-Na-anti-correlation. This was taken as sufficient justification of the 



underlying basic assumption. Very similar conclusions were obtained by Cavallo ct al. 



1998) 



Within this approach, one cannot investigate the possibility or necessity for even 



deeper m ixing, which would affect the hydrogen/helium structure of the models. Sweigart 



(1997a) has renewed the interest in such deep mixing by connecting the problem of 



ho rizontal branch morphology with that of observed anomalies, as previously suggested 
by Langer fc Hoffman (1995) . If the mixing leads to severe helium enhancement in the 
envelope, increased luminosities and stellar winds result, such that the star will populate 
the blue horizontal branch (HB), while it remains a red HB star without the additional 
mixing. 

While Denissenkov fc Weiss (1996) did not investigate the effect of helium transport, 
Sweigart (1997a) did not follow the evolution of the participating isotopes to compare 
with observations. In the present paper, we therefore attempt to close this gap by comput- 
ing full evolutionary sequences which include deep diffusive mixing and by investigating 
abundance anomalies using these self-consistent models as background models. In par- 
ticular, we want to answer the question how much helium enrichment of the envelope is 
necessary or allowed to achieve or to keep consistency with observations. 

In Sect. 2 we will discuss the nucleosynthesis aspects of the problem and review the 
observational status of the global 0-Na anticorrelation, which is a powerful tracer of 
the transport processes in the red giants. In Sect. 3 we will present and discuss the 
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Fig. 1. Observed [Na/Fe] vs . [0/Fe] for galactic globular cluster and field giants. Refer- 



ences for the data are: Fi eld: [Shetrone (1996a)|, ^hetrone (1996bl M3: [Kraft et al. (1992 



M13: iKraft et al. (1997)| , M15: [Sneden et al. (1997)t M92: [Shetrone (1996a)| , uj Cen: |Nor 



ris fc Costa (1995) . In lo Cen, black and white squares correspond to stars with [Fe/H] 
higher and lower than -1, respectively 



evolutionary models. After that, the predictions for the abundances based on our mixed 
models and post-processing nucleosynthesis will follow, before the conclusions close the 
paper. 



2. The O-Na anticorrelation: observational status and nucleosynthesis 
arguments 



Among the chemical anomalies observed in red giant atmospheres over the past two 
decades, the variations in oxygen and sodium have a special status. Indeed, the so-called 
global oxygen-sodium anticorrelation appears to be a common feature to all the globular 
clusters in a wide range of metallicity (-2.5<[Fe/Hl<-l), for which detailed abundance 



analysis of the brightest giants have been carried out (|Gratton & Ortolani 1989 


; Bneden 


2t al. 1991; 


Drake et al. 1992; 


Brown & Wallerstein 1992|; 


Kraft et al. 1993; 


Armosky et al. 


1994; 


Minniti et al. 199e : 


Pilachowski et al. 1996; 


Kraft et al. 1997; 


Kraft et al. 1998|).This 



pattern shows up in "normal" monometalhc clusters (M3, M4, M5, MIO, M13, M15, M92 , 



NGC 7002) as well as in the multi-metallicity cluster uj Cen (Paltoglou & Norris 1989; 
Norris fc Costa 1995| ). Fig. | summarizes the present observational status concerning 



the 0-Na-anticorrelation. But more importantly, this feature, and the dependence of the 



Na enhancement and O depletion on the red giant evolutionary state ( Pilachowski et al 



1996; Kraft et al. 1997) can be explained straightforwardly in the deep mixing scenario 
( Denissenkov fc Weiss 199^ . 



Regarding the Mg and Al anomalies, the situation is very different, and the observed 



Mg-Al anticorrelation (Shetrone 1996a) requires a combination of the deep mixing and 
primordial scenarios (Denissenkov et al. 1997; Denissenkov et al. 1998| ). Ad hoc assump- 
tions would be needed to obtain the observed ^^Mg depletion within the low mass Al-rich 
giants themselves ( ^hetrone 1996b ). Since a low energy resonance in the ^^Mg(p,7)^^Al 
reaction remains undetected (Angulo et al. 1999) or even can be excluded ( [PoweU 1999| ), 
"exotic" models are required to episodically increase the temperature of the hydrogen- 
burning shell up to values as high as ^ 70-85 MK (while canonical models reach a 
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maximum temperature of only 55 MK) in order to deplete Mg at the expense of ^"^Mg 
(Langer et al. 1997; Zaidins & Langer 1997; Fujimoto et al. 1999) inside the low mass 
giants as apparently being the consequence of the results by Shetrone (1996a and 1996b) 
for M13. Lately, however, Ivans et al. (1999]| found that Mg and Al abundance variations 
in M4 can be explained completely by the idea that the Al-enhancement is due to a de- 
struction of the Mg-isotopes ^^Mg and ^^Mg (cf. Fig. |^). Even in this case, a significant 
increase of the initial abundance of ^^Mg is required (Denissenkov et al. 1998). Since the 
Mg-Al ant icor relation cannot be explained by the deep mixing scenario alone, we do not 
consider it further in this paper. In fact, all explanations brought forward up to now 
being rather exotic and complicated, one should also wait for additional observational 
support for its existence (as well as for the isotopic ratios) and relation to other stellar 
properties before advocating any explanations. 



Ig L/La = 2.ai 




0,00 0.05 0,10 0.15 0.20 0,25 
6m 



Fig. 2. Abundances of some isotopes participating in the CNO-, NeNa- and MgAl-cycle as 
functions of the relative mass coordinate {5m = at the bottom of the hydrogen burning 
shell and 6m = 1 at the bottom of the convective envelope) in the red giant model from 
which detailed nucleosynthesis calculations with extra mixing start. ^^Na is emphasized 
by a thick solid line. Also shown is the hydrogen abundance profile, normalized, however, 
to 10"'^ to accomodate it on the same scale 



Last but not least, the morphology of the global 0-Na anticorrelation (Fig. ||) bears 
crucial clues on the mixing process. First, its extension to very low oxygen abundances 
([0/Fe] < —0.45) is entirely due to the contribution of M13 red giants. It is this "second 
para meter" globular cluster that has the fastest rotating blue horizontal branch (HB) 
stars. Peterson et al. (1995) found six stars having v sini > 30 kms^^. In the same paper 
red HB stars in M3 (the cluster forming a classical "second parameter" -pair with M13) 
have been found to possess smaller projected rotational velocities, from 2 to 20 kms^^. 
This indicates a relation between rotation and mixing. 

Concerning this extension along the horizontal axis, our experience in modelling the 
global anticorrelation shows that it depends primarily on the mixing rate Dmix or, more 
precisely, on the product "mixing rate x mixing time" . We will come back to these 
arguments in Sect. 4, where our nucleosynthesis predictions in deeply mixed stars will be 
compared to the observations of [0/Fe] versus [Na/Fe]. 
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Secondly, the extension of the global anticorrelation along the vertical axis (—0.4 < 
[Na/Fe] < 0.6 in all clusters except ui Cen) tells us mostly about the depth of the 
additional mixing. In Denissenkov & Weiss ( |1996 , Fig. 4) and in Denissenkov et al 



(1998, Fig. 1) it has been shown that on approaching the hydrogen burning shell the 



Na abundance displays two successive rises (see Fig. |). The first rise results from the 
reaction ^^Ne(p,7)^'^Na, whereas the deeper one is produced in the NeNa-cycle by the 
partial consumption of ^"^Ne, which is much more abundant than both ^^Ne and ^^Na. The 
size of the vertical extension of the global anticorrelation implies that additional mixing, 
whatever it is, does not penetrate the second step (rise) in the Na abundance profile 
where H starts to decrease. Otherw ise, observed [Na/Fe] values would be much larger 
than they areQ. [Cavallo et al. (1998) have investigated the dependence of the abundance 
profiles on metallicity (and mass). Their results confirm our arguments completely. Only 
for stars of near-solar metallicity very deep mixing with significant heliu m enrichment 
but without st rong sodium enhancement could be possible (see also Fig. 5 of Denissenkov 
fc Weiss 1996| ). However, the clusters under discussi on (e.g. M13) are metal-poor. Le t 
us note in Fig. |l] the different behavior of field stars ( Bhetronc 1996a , shetrone 1996b ). 
In this population, the 0-Na-anticorrelation is not present ( Gratton et al. 2000 ). This 
indicates that the deep mixing does not penetrate the region where ON-burning occurs, 
and may reveal possible environmental effects on its efficiency. 

To prepare Fig. ^ we applied a nucleosynthesis code to a red giant model with surface 
luminosity log L/Lq = 2.21 from which our nucleosynthesis with additional deep mix- 
ing calculations started (Sect. 3). The considerable growth of the ^^Al abundance with 
depth is due to our "non-standard" assumption of an enhanced initial abundance of the 
^^Mg isotope ([^^Mg/Fe] = 1.1) and of the thousandfold enhanced rate of the reaction 
^^Al^(p, 7)^^Si (for details see Denissenkov et al. 1998). These ad hoc modifications were 
needed to explain the observed Al enhancements (Denissenkov et al. 1998), but have no 
influence on the results of the present paper. From Fig. ^ we see immediately that if 
extra mixing penetrates down to layers, say, at Sm w 0.06, this will result in an enrich- 
ment of the red giant's envelope in N, Na and Al and in its impoverishment in C, O, and 
^^Mg. Again, ^^Mg remains unchanged due to the relatively low temperatures reached 
in such a star. These results were recently confirmed by Palacios et al. (1999) in models 
using the reaction rates recommended by NACRE ( Angulo et al. 19991 ). 



3. Red giant evolution with deep mixing 

To produce background models for the nucleosynthesis post-processing we have evolved 
stellar models under the assumption of additional deep mixing after the RGB bump. All 
sequences were started at the same initial model, which consisted of an 0.8 Mq star of 
initial composition Y — 0.25 and Z — 0.0003, and which had been evolved (canonically) 
up to the luminosity of the bump, i.e., logL/L© = 2.21 (at this point, its mass is 
0.798 Mq). The envelope helium content has increased to 0.256 (in mass fraction) due to 
the first dredge-up. The input ph ysics of the Garching ste llar evolution code (used here) 



is up-to-date (for a summary, see Denissenkov et al. 1998 ), but atomic diffusion has not 
been included in the computations. 

The additional deep mixing between the convective envelope and some point inside 
the hydrogen shell has been implemented in the same general line as in our previous 
papers on this subject, that is, as a diffusive process with parameterized values for the 
diffusive constant I?mix, which is the same for all elements, and for the penetration 
depth. The values used for Dmix are guided by the results of our earlier papers, and 
agree with estimates based on rotationally induced mixing theories. We refer the reader 



to Denissenkov & Weiss (1996) for details of this approach. We use the normalized mass 



^ Actually some metal-rich giants in lj Cen show extremely high [Na/Fe] (up to 1 dex), prob- 
ably indicating a penetration of the mixing to the second rise of Na. This cluster is the only one 
where some giants exhibit surface enrichm ent of Na produced f rom both ^"Ne and ^^Ne. It also 
has one of the bluest horizontal branches ( Whitney et al. 199'j ). 
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coordinate 6m introduced therein, which is at the bottom of the hydrogen shell (usually, 
where X = 10^^) and 1 at the bottom of the convective envelope. The choice of this mass 
coordinate allows accurate interpolation between a small number of background models in 
the nucleosynthesis calculations (see Cavallo et al. 1998 for a similar approach). The shell, 
in this coordinate, is located below Sm w 0.10. The depth, down to which the diffusive 
mixing should occur, we denote as (Sm„iix- Obviously, due to the lack of solid theories, 
one could also choose, for example, purely geometrical scales to define the penetration 
depth (Boothroyd & Sackmann 1999b). Contrary to our earlier papers, the criterion for 
penetration is not determined by a fixed value for Sm„iix chosen before the calculations, 
but is related to the decrease in hydrogen content within the shell (relative to the surface 
or convective envelope abundance Xcm), expressed as a free parameter AX. We have 
investigated several different prescriptions for the penetration criteria and found a great 
sensitivity of the mixing on these prescriptions, which are 

1. find that Sm^i^ in the initial model, where X — Xc^v — AX, and mix to the same 
Jrrimix in all subsequent models; 

2. as 1., but 6m = is defined as the point where X = Xonv/2 (instead oi X = 10^"*); 

3. as 1., but the diffusion constant Dmix is decreasing exponentially from the maximum 
value Do for 6m > 0.10 to i^mix ~ 5 • 10~^I?o at Sm^i^ 

4. always mix to the point, where X — Xonv — AX 

Except for method 3, these schemes were chosen to be as simple as possible and to be 



similar to the one by Sweigart (1997a) 




0.10 



Fig. 3. Hydrogen profiles in the normalized 6m coordinate in models with deep mixing 
during the evolution along the RGB. Mixing parameters were Dmix = 10^ cm^ s~^; AX = 
0.20, method 1. The vertical line is the mixing depth ^m,nix as defined in the first model. 
See text for more explanations 



Keeping the relative mass coordinate fixed, up to which mixing should occur, im- 
plies that changes in the hydrogen profile in the shell (usually steepening in the course 
of evolution) or in the extend of the convective envelope influence the mixing. As an 
illustration of the "movement" of the profile in this coordinate we display in Fig. |^ an 
example (i^mix — lO^cm^s"^; AX = 0.20; method 1). The smooth solid line is the 
initial model. The next model (at log L/Lq — 2.288) is the left-most line. From there, 
the model profiles shift to the right again. The solid line with clearly reduced hydro- 
gen abundance throughout the envelope corresponds to a model close to the RGB-tip 
(logL/L© = 3.539). In this phase the burning time at the bottom of the mixed region 
becomes short enough to lead to "bottom-burning" of the envelope. An extended RGB- 
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phase with luminosities drasticaUy increased above the canonical RGB-tip luminosity of 
l og -Ltip 7-^0 = 3.33 is the result; the final value being logLtip/-^© — 3.8. Since mass loss 



(Reimers f975, with 77 = 0.3) was taken into account, the total mass at the tip decreases 
in such cases below 0.6 Mq. ^m-mix is indicated by the vertical line and is 0.047. 

Using method 2 instead, Sm = is defined as the point where the abundance of 
hydrogen has dropped to half the surface value of the same model (Smmix = 0.0066 in 
this case). This prevents, obviously, any penetration to shell regions with lower hydrogen 
content. As in the previous case, the bottom of the envelope is burnt at the end of 
the RGB evolution. However, as soon as the hydrogen abundance approaches Xonv/2 
(guaranteed as long as mixing is not quasi- instantaneous) , the mixing criterion inhibits 
further mixing. For this reason X remains constant for Sm < 0; in the final models 
^cnv — 0.6402 and the hydrogen profile always has a finite step in the shell. In this case, 
the luminosity rises to logLup/ Lq — 3.53. 

The steps visible in the chemical profiles of Fig. || are due to the diffusion criterion 
applied to the numerical grid, because no interpolation to the exact value of (5mmix had 
been done. We verified that the results do not depend on the grid resolution, which 
was increased by a factor of 10 in the shell in part of the calculations. Only the steps 
got smaller and more numerous. To avoid such steps, we introduced a varying diffusion 
constant (method 3) motivated by recent results of Denissenkov fc Tout (2000)] . They 



have proposed a physical mechanism for extra mixing in red giants which quantitatively 
interprets all the known star-to-star abundance variations in globular clusters. This is 
Zahn's mechanism ( ^ahn 1992 ; Maeder & Zahn 1998) which considers extra mixing in a 



radiative zone of a rotating star as a result of the joint operation of meridional circulation 
and turbulent diffusion. This process was already advocated by [Charbonnel (199*5) 



Do exp 



CD 



< Sm < Smo (1) 



to 

explain the low carbon isotopic ratios and lithium abundances in field Population II giants 
and to lower the ^He yields by low mass stars. Denissenkov & Tout report that the mixing 
rate does not vanish abruptly at a particular depth but instead it dies out gradually on a 
rather short depth range approximately between Sm = 0.10 and Sm — 0.06 ~ 0.07. This 
explains the following choice of an exponential decline approach for -Dmix: 

-Dmix = Dq- Sm> Smo 

Smo — Sm 

(^TOinix - Smo^ 

where Smo — 0.10 was used for the beginning of the decline and (5m,nix is the mixing depth 
coordinate as defined in method 1. Using cjj = 10 ensures that £'mix('5TOmix) ~ 5-10^'^£'o- 
The resulting evolution (Fig. ^) is similar to that of the case shown in Fig. ^, but the 
profiles are smooth; mixing parameters are identical. 

We finally note that method 4, applied straightforwardly, leads to a complete burning 
of the entire envelope for Dmix > 5 • 10^ and/or AX > 0.10. However, this we consider 
to be an artefact, which is easy to understand: Since due to the mixing the hydrogen 
abundance in the envelope is reduced, the critical point down to which mixing should 
occur, is moving inwards. Therefore layers of even lower hydrogen content are mixed with 
the envelope, and the critical point moves to even deeper regions. Only at the point where 
the burning time-scale is shorter than the mixing time-scale and therefore the surface 
hydrogen content no longer is able to adjust to the burning, this process is stopped. A 
way out of this situation is to ensure that the mixing does not lead to a sharp step in 
the shell's hydrogen profile. This way, the critic al point can be k ept within the mixed 
layers. While the mixing procedure described in Sweigart (1997a)| appears to follow the 



straightforward approach, Sweigart (private communication, 1999) in fact used a more 
complicated method to keep the hydrogen profile even in the presence of mixing. Our 
method 3 qualitatively has the same effect. 

To summarize, method 4 appears to be unphysical; method 3 is based on the physical 



picture by Denissenkov & Tout (2000) , but makes additional parameters necessary. Meth- 
ods 1 and 2 are the most straightforward choices leading to composition profiles similar 
to method 3; method 2 differs from 1 in that it avoids complete mixing and burning of 
the envelope during the very last phases of RGB evolution. 
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0.10 



Fig. 4. As Fig. ^, but for an exponentially declining diffusion coefficient Dmix (Eq. 
The results shown arc from case D in Tab. 1 

Table 1. Parameters of calculations presented in Figs. and final values of the stel- 
lar mass {M{/Mq) and luminosity (logLf/LQ), and of the helium mass fraction in the 
envelope (ycnv)- The final luminosity Lf is always very close or identical to the RGB-tip 
luminosity Ltip. The mass loss in sequences B and C has been reduced by a factor 20 
compared to the other ones. The penetration depths, expressed in the normalized mass 
coordinate 6m corresponding to the three AX values are 0.060 (0.05), 0.054 (0.10), and 
0.047 (0.20). Anix is in units of cm^ s~i 



case 




AX 


method 


M{ 


Y 

cnv 


log if 


S 


0.0 






0.689 


0.256 


3.31 


A 


5 ■ 10* 


0.05 


3 


0.586 


0.270 


3.79 


B 


5 ■ 10* 


0.10 


2 


0.792 


0.284 


3.30 


C 


10^ 


0.20 


2 


0.567 


0.360 


3.53 


C 


10^ 


0.20 


2 


0.771 


0.360 


3.56 


D 


10^ 


0.20 


3 


0.562 


0.351 


3.77 



Obviously the details of the mixing procedure influence the resulting evolution to quite 
a significant extent. Since at present we are far from providing a solid physical approach 
(which could allow, for example, for a diffusion speed varying both in space and time), 
we cannot predict the true evolution of a star experiencing deep mixing. However, the 
calculations are needed only to provide background models with varying degrees of helium 
mixing into the envelope, and it is of no importance how this is achieved in detail. We have 
calculated 25 different sequences, varying method and parameters. The cases selected for 
Tab. |ll are representative for the range of results we obtained, which are summarized in 
Figs. ^, ^, and 0. Case D of Tab. is the one also shown in Fig. ^. We add that the 
amount of mass loss has no influence on the mixing properties. Neither does a gradual 
switching-on of the additional mixing during the first few models. We have performed 
some comparison calculations with a completely different code (the Toulouse-Geneva 
code; Charbonnel et al. 1992| ). While the results differ in details, the gross properties 
are the same. The differences we ascribe to details in the mixing procedures and the 
implementation of diffusion. 

The effects on the evolution, displayed in Figs. || ~ 0, are qualitatively as expected 
from the work by Sweigart (1997a). The increase in the surface helium content is quite 
dr amatic in cases C and D (fast and very deep mixing). However, it is not as large as 
in Sweigart (1997a) , shown, for example, in his Fig. 2, where for AX — 0.20 a value of 
Y ^ 0.42 was reached. Also, in contrast to Sweigart's result, in all our calculations the 
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helium enrichment of the outer envelope tends to level off with progressing evolution. 
This might be ascribed again to differences in the mixing scheme details, as we find these 
differences also in the post-processing models presented in the next section (cf. Figs. ^ 
and We also find that the luminosities can get extremely high (cases A and D) with 
high mass loss as the consequence and a beginning turn-away from the RGB before the 
He-flash sets in. The beginning of such an evolution might be recognised, too, in Fig. 3 
of sweigart (1997b) in the case of deepest mixing. We also note that mixing method 2 
(cases B and C) results in loops in the HRD, which depend on the occurrence of mixing 
episodes. This is, for example, visible in the non-monotonic luminosity evolution in Fig. ^. 



4.0 




3.70 3.68 3.66 3.64 3.62 3.60 

log Trf, 

Fig. 5. HRD of four mixed sequences of Tab. Q and the unmixed canonical one (S) for 
comparison. The linotypes refer to the cases listed in Tab. |l| and are given in the top- 
left corner. The large plot shows the upper RGB evolution only, the inset the complete 
evolution from the ZAMS on. 



0.36 



0.34 



0.32 - 



0.30 - 



0.28 



0.26 



0.24 




4.0 



3.0 

log L/U 

Fig. 6. Surface helium abundance (in mass fraction) as a function of luminosity for the 
same sequences as in Fig. || A coarseness of 0.001 is due to the limited number of digits 
in the output 
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1 .290 



1 .292 



1 .291 
age ( 1 0'°yr) 

Fig. 7. Surface helium abundance as a function of age for the same sequences as in Fig. || 



The hfctime on the RGB is in aU cases prolonged (Fig. |^) by 1 or 2 Myr, in some 
cases up to 4 Myr (5-10% of the lifetime after the bump). The evolutionary speed is 
influenced mainly immediately after the onset of the additional mixing (after the bump) . 
From the initial model to one at log L/Lq « 2.5, the time increases from 7 Myr (case S) 
to about 10 Myr (A & B) and 13 Myr (C & D). Thereafter, it is becoming comparable 
again in all cases. Such an effect could possibly be seen in luminosity functions, but it 



should be most prominent only in a limited luminosity range. VandenBerg et al. (1998) 
have argued that the luminosity function of M30 could be evidence for rapidly rotating 
cores of giants. 

In some calculations, in particular those employing the exponentially declining diffu- 
sive speed (method 3), an interesting effect appeared. Although the mixing for moderate 
penetration depth and mixing speed remains small, the luminosity of the models rises 
to extreme values. As an example, in case A {AX = 0.05, Z^mix — 5 • lO^cm^s"^), 
Ycnv = 0.270 (an enrichment of only +0.014) and logL/L© = 3.79 were reached at the 
tip of the RGB (see also Fig. H). An inspection of the models shows that part of the ex- 
tended region of almost homogeneous composition, which is achieved due to the effect of 
the additional diffusion, becomes hot enough for significant hydrogen burning. In Fig. ^ 
we display the H-profiles of selected models in this phase. The first one (top line) is at 
\ogL/LQ — 2.73 (Vcnv — 0.265). In this model, the point were the energy production due 
to hydrogen burning exceeds 10'^ ergg~^ sec~^ for the first time, is very close to the com- 
position step. In the next model (logL/L© — 3.12; Y^nv — 0.269) this point has shifted 
to 6m — 0.274 and is constantly progressing outward until it reaches 6m — 0.80 in the 
most advanced models at the bottom of the figure (logL/L© = 3.68; ycnv = 0.270). This 
burning of the plateau constitutes a broadening of the hydrogen shell and delivers extra 
luminosities. In fact, 50% of the total luminosity of the models with the plateau value 
around X = 0.40 are generated at 6m > 0.06, that is outside the inner composition 
step. Therefore, the luminosity in excess of that of an ordinary star at the tip of the 
RGB (around log L/Lq = 3.30) can completely be ascribed to the plateau burning (note 
that diffusion cannot keep the burning plateau homogeneous with the outer regions). 
The effect we observe here is probably due to our mixing description, which sets Dmix to 
the maximum value outside 6m = 0.10. If our mixing description is realistic, this would 
imply that one could get very high luminosities at the RGB-tip without extreme helium 
mixing. The time spent at luminosities above the standard TRGB brightness is only 
10^ yrs and therefore observation of such a superluminous star is rather unlikely. Due 
to the extreme overluminosities, the Reimers mass loss formula leads to stellar winds of 
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order lO-'^Me/yr and a final mass of 0.58 Mq (Mc = 0.517 Mq). After the helium flash, 
this star will populate the blue part of the horizontal branch. 



0.8 




0.0 iL ^ ^ \ ^ . ^ \ ^ ^ ^ \ ^ . ^ \ ^ ^ ^ 

0.0 0.2 0.4 0.6 0.8 1.0 

5m 



Fig. 8. Hydrogen profiles of selected models at the tip of sequence A (AX = 0.05, 

Dmix = 5-10* cm^ s""'^, mixing method 3). The point where energy generation by hydrogen 
burning exceeds 10^ erg/ (gs) first, is indicated in several models by symbols. 
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Fig. 9. The global anticorrelation and three theoretical curves (labelled by the value 
of r'mix used in the corresponding calculation); all three curves were calculated with 
the penetration criterion 1 and AX = 0.2; unmixed red giant models were used as 
background. Observational data (omitting oj Cen) are as in Fig. 



4. Nucleosynthesis in deeply mixed stars 

The detailed nucleosynthesis calculations we present now have been performed in a post- 



processing way as in our previous works (Denissenkov & Weiss 1996; Denissenkov et al 
1998| ). From an evolutionary sequence three red giant models were selected. The starting 



one was the same one as for the full evolutionary calculations discussed in the previous 
section, i.e. a model at the bump (logL/L0 w 2.2) in which the hydrogen burning shell 
had recently crossed the H-He discontinuity left by the base of the convective envelope on 
the first dredge-up phase. The final one was a model near the RGB tip (logL/L© « 3.3), 
the second one having a luminosity intermediate to those of the starting and finishing 
models. Distributions of T, p and r with Sm in these three "background" models were used 
for interpolations in log L during the nucleosynthesis calculations. Further details about 
our post-processing procedure can be found in Denissenkov fc Weiss (1996)| The network 



of nuclear kinetics equations was the smallest one of those considered in Denissenkov et al 



(1998). It takes into account 26 particles coupled by 30 nuclear reactions from the pp- 
chains, CNO-, NeNa- and MgAl-cycle. The additional mixing is modelled by diffusion 
with a constant coefficient -Dmix- We recall that we allow for mixing prescriptions and 
parameters in these calculations different from those for which the background models 
have been obtained. 

For the comparison with observations we have preferred the "global anticorrelation" 
of [0/Fe] versus [Na/Fe] for the reasons detailed in Sect. 2. In Fig. ^ and Fig. ^ it is 
plotted for globular clusters according to the latest observational data. 

In Fig. ^ theoretical dependences of [Na/Fe] on [0/Fe] obtained in the post-processing 
way for three values of the diffusion coefficient are shown. This first set of calculations 
was performed with unmixed background red giant models like in our previous papers 
but with mass loss taken into account. The depth of additional mixing was determined 
according to the penetration criterion 1 (Sect. 2) with AX — 0.2. The mass loss rate 
M was estimated wit h pR,eimers (1975) formula in which the parameter value 77 = 0.3 
was adopted. In Fig. |10| the resulting envelope He abundances are shown as functions 
of log L/Lq. The mixing depth in the starting model which was kept constant during 
the nucleosynthesis calculations was Sm^i^ — 0.047. Such a value of (Jm^ix allows some 
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Fig. 11. As Fig. but for calculations using mixed red giant models as background. 
The solid curves were calculated with the penetration criterion 1 and AX ~ 0.2 and the 
dashed ones with constant mixing depth (5m„iix = 0.04 



(modest) penetration of the second Na step (see Sect. 2) by the mixing which results in 
an upward steepening of the theoretical dependences of [Na/Fe] on [0/Fe] by the end of 
the RGB evolution (Fig. The maximum He enrichment achieved in the envelope in 
this set of calculations is Al^nv ~ 0.10 (Fig. [l0| ). 

In the next nucleosynthesis calculations, the results of which are presented in Figs. ^ 
and we have used mixed background models from sequence C of Tab. |l|; for this 
sequence the same mixing parameters as in case C (Tab. ^ but a reduced mass loss rate 
has been used: the parameter 77 was divided by 20 in order to take into account a reduction 
of M for low Z (Maeder 1992). The reduced mass loss rate affects only the final mass, 
but not the helium enrichment. These models were evolved along the RGB with effects 
of the He mixing on stellar structure parameter distributions fully taken into account in 
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the stellar evolution code (Sect. 3). Sequence C was chosen because out of the sample 
cases listed in Tab. |l| it has the highest degree of helium enrichment, in contrast to the 
previous, unmixed, background models. In this second set of nucleosynthesis calculations 
we repeat the case of the first set (Fig. i.e. mixing down to AX = 0.2 (solid lines), but 
also add two computer runs (dashed lines) in which the mixing was chosen to be so deep 
that a high enrichment in He of the envelope was guaranteed. We label these calculations 
by (5minix = 0.04, which is the penetration depth needed to mix down to AX — 0.37. 

Comparison of the results obtained in the two sets of calculations allows to draw the 
following conclusions: 

— mass loss is practically unimportant for this study (at least within the prescriptions 
and variations used in the various calculations); 

~ making use of mixed background models instead of unmixed ones does not seriously 
affect the theoretical dependences of [Na/Fe] on [0/Fe] (compare the solid curves in 
Figs. ^ and |ll|); therefore the details of the mixing prescription used for the back- 
ground models are not significant for the nucleosynthesis results. 

— the total He enrichment of the convective envelope calculated in the post-processing 
way agrees very well with the final envelope He abundance obtained in the full evo- 
lutionary calculations with additional deep mixing; 

— values Aycnv > 0.15 were obtained only in the two computer runs with the depth 
Smmix — 0.04, but in these cases additional mixing penetrated so deeply that it 
resulted in the [Na/Fe] on [0/Fe] dependences evidently inconsistent with the obser- 
vations (dashed curves in Figs. ^ and |lj). 

A simple inspection of Figs. p[]l^ allows the conclusion that the global anticorrelation 
of [0/Fe] vs. [Na/Fe] as a whole and especially the [Na/Fe] values in its low oxygen abun- 
dance tail certainly rule out any hypothesis about an increase of more than AFenv ~ 0.10 
in the envelope He abundance of globular- cluster red giants. A physical reason for this 
constraint is the above-mentioned inability of additional mixing to penetrate the second 
Na abundance rise lying at Sm < 0.06 0.07 as hinted by the observed global anticorre- 
lation. In the starting models Sm > 0.07 corresponds to AX < 0.05 and, therefore, the 
envelope He enrichment is not expected to be much larger than AY^nv ~ 0.05. 
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Fig. 13. [Na/Fe] and [0/Fe] for giants in globular clusters. The symbols identify indi- 
vidual clusters and have the same meaning as in Fig. ^ Overlaid are the five theoretical 
predictions shown in Fig. |l^, labeled by the values of -Dmix- In the lower panel the two 
higher values relate to two lines each. 



5. Discussion 



If the 0-Na-anticorrelation observed in many globular cluster red giants is indeed due 
to a deep mixing process beyond the standard effects taken into account in canonical 
stellar evolution theory, the question is justified whether this deep mixing might affect 
the H-He-profile as well. In this case, consequences for the red giant evolution including 
phases of enhanced luminosities and mass loss could result. We have, therefore, discussed 
both the evolutionary and nucleosynthetic effects quantitatively by performing stellar 
evolution calculations including deep mixing and post-processing nucleosynthesis models 



(the latter as we did in our earlier papers Denissenkov & Weiss 1996; Denissenkov et al. 



1998). 



From arguments depending only on nucleosynthesis we could already infer that for 
temperature profiles typical of hydrogen-burning shells mixing of appreciable amounts 
of helium can only be achieved if the second Na rise is penetrated. This, however, leads 
to oxygen and sodium anomalies exceeding those observed (with the exception of a few 
stars in uj Cen, a multi-metallicity, untypical cluster). In terms of our normalized mass 
coordinate (defined such that 5ra = at X = 10"'* at the bottom of the shell) this puts an 
observationally constrained limit for the maximum mixing depth of ^m„iix > 0.06-^0.07. 
Our complete models confirm this argument: helium enrichment in excess of Al^nv ~ 0.05 
due to deep mixing can be ruled out for those stars with Na-O-anomalies as observed in 
clusters such as M15, M92, M3, and even M13 which presents one of the most extended 
blue horizontal branch. 

The details of the surface abundance history along the RGB depend on the details 
of the deep mixing process and thus on its nature which we did not attempt to specify 
here. In particular, the field-to-cluster differences must be understood (especially the fact 
that field giants do not present the 0-Na anticorrelation ( Gratton et al. 2000| ) , indicating 
a deeper and more efficient mixing in their globular cluster counterparts) and seem to 
point out a non-negligible impact of environmental effects on the extra mixing efficiency. 
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However, we can compare the observations with the histories of Na and O abundance 
anomahes predicted by our simple mixing prescriptions in order to get constraints for 
a sohd physical model. This is done in Fig. which shows surface abundances for 
several clusters as a function of brightness, i.e. progressing evolution. Also shown is 
the theoretical prediction of the five calculations displayed in Fig. Although the 
observational data are very few (the uncertainty in the abundances is of order 0.2 dex), 
and information for stars before or at the bump is available only for M13, some effects 
can be recognized, nevertheless. 

All clusters show the whole spread of Na abundances between the canonical case 
without extra mixing (they would lie on a horizontal line) and that as obtained from 
the calculations with less deep mixing (solid lines). As already mentioned, very low O 
abundances are only found in M13 giants, and the observed spread for this element is a 
signature of the mixing rate. For both elements the abundance anomalies are limited to 
values predicted by models with extra mixing not penetrating the second ^'^Na rise in the 
hydrogen shell. Stars with intermediate anomalies we interpret as being due to mixing 
penetrating less deeply into the hydrogen shell. They could be reproduced with properly 
adjusted mixing parameters. 

Over the small brightness interval for which we have data (for all the clusters except 
M13 only the brightest giants are accessible), no significant abundance evolution is rec- 
ognizable. The increase in [Na/Fe] at the RGB tip obtained in all calculations labeled 
AX = 0.20 is not visible in the observations, which may be taken as indication that the 
physical reasons (e.g. rotation) for the additional mixing have lost their importance or 
even vanished. Such a possible time dependence of the extra mixing has not been taken 
into account in our simple mixing prescriptions, but should result from more physically 
motivated models. (We recall that in Denissenkov fc Weiss (1996)| the carbon evolu- 



tion could be reproduced with constant mixing parameters, however.) We will therefore. 



in a forthcoming paper, use the model by Denissenkov & Tout (2000), which includes 



temporal changes in the diffusion constant due to angular momentum transport. 

In our models, the largest abundance changes take place at the onset of the additional 
mixing, i.e., early on the RGB, after the bump (around My ~ —0.5). This might be 
visible in the M13 data. The low brightness group has normal abundances, but above 
My ~ —0.4 strong anomalies already appear, reminiscent of the steep increase in the 
calculations mixing deeper (dashed lines). In the case of Na no further enhancement 
is visible (the range of values does not varies along the RGB), while O seems to get 
depleted further. In spite of the incomplete data at hand, the abundance anomalies seem 
to develop rather early and within a narrow brightness range, but then do not increase 
any more. This general behaviour is consistent with our theoretical predictions - though 
not exactly reproduced ~ and again points to not too deep mixing at moderate speed. The 
calculations with extreme helium mixing would predict the largest O depletion and Na 
enhancement all along the RGB; the absence of such stars can therefore not be explained 
with a selection effect working against the shortest lived stars at the tip of the RGB. 
The fact that M13 seems to show anomalies already before the bump has to be taken 
with care, because we compare here only with one stellar model which has a metallicity 
almost a factor of 10 smaller than that appropriate for M13. At that metallicity, the 
bump would occur about 1 mag earlier (and our initial model is about 0.2 mag brighter 
than the end of the bump). We stress the fact that the present observational data do not 
allow a detailed comparison with the abundance evolution, but only rough qualitative 
statements. 

Although Fig. |l^ again demonstrates how the abundance evolution depends on the 
assumption about the additional mixing process, our calculations also show that the 
nucleosynthesis argument given above is valid for all mixing descriptions tested and 
therefore model-independent. From this we predict that red giants exhibiting the observed 
Na-O-anomaUes do not have envelopes enriched in heUum by much more than AY^nv ~ 
0.05, which is comparable to the general uncertainty in our knowledge about the helium 
abundances in such stars. 
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Concerning the brightness reached during the RGB evolution, we showed that for 
specific mixing prescriptions large excess luminosities can be achieved at the end of the 
evolution without simultaneous mixing of large amounts of helium. As the reason we 
could identify the burning of the outer parts of the hydrogen shell where some H-He- 
mixing had happened in earlier phases. These regions become hot enough for significant 
hydrogen burning on such short time-scales that the deep mixing process is not able to 
mix the products any longer into the convective envelope. Such models experience very 
strong mass loss under the assumption of the continuing validity of a Reimers-type stellar 
wind and finish the RGB phase with total masses below 0.6 M© and envelope masses of 
only 10% of this. They could be candidates for blue HB stars and would link observed 
ab undance anomalies, deep m ixin g and the second parameter problem, as suggested 
by Langer fc Hoffman (1995) and Sweigart (1997a)| . They would avoid the problem of 



overproducing 0- Na-anomalies as would result from the helium mixing investigated by 



Sweigart (1997a) 



However, we repeat our warning that the details of the evolutionary consequences of 
deep mixing depend crucially on the mixing process and history. Mixing speed and depth 
are both important for the amount of helium mixed and for the result of the competing 
processes "mixing" and "burning" and their time-scales. The mixing depth is also linked 
to the criterion for deep mixing, for which we have investigated several simple recipes. 

To conclude, we need a solid physical picture for the deep mixing process in order 
to be able to investigate its effect on red giant evolution further. Presently, we can only 
point out some interesting possibilities - such as the overluminosities - and derive model- 
independent features, such as our main conclusion that the observed anomalies of oxygen 
and sodium rule out strong helium enhancement and therefore very deep mixing. 
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